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Abstract 

We analyze the electromagnetic pion and kaon form factor by including radiative and higher-twist 
effects within the framework of resummed pQCD in the space-like region. We focus on the transition 
from the perturbative to non-perturbative behavior in the phenomenological intermediate energy 
regime. Using a modified "A;t" factorization scheme with transverse degrees of freedom, we evaluate 
the non-perturbative soft contributions as distinct from the hard contributions, ensuring no double 
counting via the Ward identity at Q 2 = 0. The soft contributions are obtained via local quark- 
hadron duality while the hard contributions rest on the well known collinear factorization theorem 
using model wave functions with modified Brodsky-Huang-Lepage type ansatz and distribution 
amplitudes derived from light-cone QCD sum rules. Our analysis shows that the perturbative hard 
part prevails for large Q 2 beyond 50-100 GeV 2 , while for low and moderate momentum transfers 
below 10-16 GeV 2 the soft contributions dominate over the hard part. Thus, we demonstrate the 
importance of including the soft contributions for explaining the experimental form-factor data. 
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I. INTRODUCTION 



During the past decade, QCD-oriented studies have been shifting steadily toward exclu- 
sive channels. For a long time the electromagnetic structure of pions has been subjected to 
numerous experimental and theoretical investigations through the study of electroproduction 
reactions. Extensive experimental studies of pion electroproduction reactions like ep^e7i + n 
or en^en^p have been carried out in the past at CERN, Cornell, DESY and more recently 
at JLab Facility . Since the mid-90s, kaon electroproduction reactions 

like A(j, K)YB and A(e, e'K)YB (A is the target, Y the produced hyperon and B the 
recoil) have also attracted renewed interest in nuclear physics at both experimental 0, [h]] 



and theoretical 



11 



121 ] level. The main ingredients for the description of electromagneti- 



cally induced kaon production are embedded in the so-called Chew, Goldberger, Low and 
Nambu (CGLN) scattering amplitudes. In case of the longitudinal component of the elec- 
tron induced unpolarized differential cross-section, the t-channel diagram dominates and (in 
certain kinematic conditions) can be factorized [si, Q as <jl = k- J-'iQ 2 ) Q(W)H(t), where k 
is a kinematic factor, T, Q, and 7i are functions of the 4-momentum transfer squared Q 2 of 
the virtual photon, the invariant mass W and the Mandelstam variable t, respectively. The 
function T{Q 2 ) implicitly contains the information about the electromagnetic form factor of 
the kaon. Note that ^{Q 2 ) is not the actual form factor, but rather a complicated function 
from which the form factor can be extracted using e.g., Chew-Low extrapolation and decon- 
volution algorithms. A precise knowledge of the form factor is of fundamental importance 
for a realistic and accurate description of exclusive reaction mechanisms and plays a key 
role in understanding the interplay between perturbative and non-perturbative physics at 
intermediate energies. Moreover, the study of form factors provides direct insight into the 
electromagnetic structures and charge distributions of hadrons as they couple with photons. 

To date, the electromagnetic kaon form factor is very poorly known and only measured 
at very low Q 2 (below 0.2 GeV 2 ) HQ. 

The status for the (quasi-free) Lambda (A) 
and Sigma (S) hyperons is even worse, i.e., there are simply no available experimental 
data. Basic quantities like the strong coupling constants Qkkn and Qky.n derived from 
purely hadronic processes or theoretical considerations are not well established and must 
be considered adjustable. Recently, however, there appeared quite large and precise data 
sets on photo-production of kaons from the SAPHIR (ELSA) [15], CLAS (CEBAF) [16] and 
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LEPS (SPring8) 17j collaborations. There is also new data on electroproduction of positive 
kaons from experiment E98-108 at CEBAF which are being analyzed at the moment. 

Keeping in mind the increasing accuracy of experimental data, an accurate theoreti- 
cal description of the electromagnetic form factors of pseudo-scalar charged mesons at in- 
termediate energies is of primal importance. To our knowledge, especially for the kaon 



there are very few theoretical works in this direction 



18 



19 



20 



2l|. In this paper, 



in addition to the pion form factor, we analyze the kaon form factor for a broad range 



of space-like momentum trans: 



light-cone QCD formalism 



22 



er. Our framework is based on resummed perturbative 



231 ] . unlike conventional approaches like "asymptotic" and 



lattice QCD or from sum rules that rely on many unchecked hypotheses. The exper- 



imental results could then be used to extract 
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ad, 
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32 



33 



34 



35 



;he various distribution amplitudes (DAs) 
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37 



38|, 



39 



40 



4l| used in the above for- 



mulation. Of course, only a handful of experimental hadron electroproduction data points 
are presently available to make definitive statements on the validity of different theoretical 
approaches. Since in almost all cases, the corresponding data points are merely concentrated 
in the very low energy region (Q 2 < 1 GeV 2 ), perturbative QCD (pQCD) has limited predic- 
tive power due to the rapidly growing magnitude of the strong coupling, as Q 2 tends to zero. 



Despite the existing plethora of literature on the predictions on electromagnetic meson 



factors based on various approaches (see e.g., Refs 



47 



48 



49 



5(1 



51 



52 



53, 



54. 



for the pion and Refs. [18 
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2l| for kaon form factors, to give a highly incomplete list of 
references), till date there is considerable amount of debate as to their exact behavior in the 
phenomenological low and moderate energies between Q 2 « 4-50 GeV 2 . Nevertheless, we 
try to give our assessment to the existing scenario and try to explain the experimental data 
first for the pion form factor, where statistics are far more decent as compared to that of the 
kaon. Then we extend our analysis to the kaon form factor, where experimental data is still 
too limited for any meaningful comparison. Hopefully, with the planned 12 GeV upgrade 
proposal of the CEBAF experiment (at JLab) in the near future, studies of intermediate 
energy QCD can prove to be fruitful. 

The standard "asymptotic" QCD is known to make successful predictions of many phe- 
nomena like dimensional scalin g, h elicities, color transparency, etc. for exclusive processes, 



as Q 2 tends to infinity [47, 



70 



71 



73]. The approach relies on the so-called collinear 
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factorization theorem 7l|, ]72|] which provides an outstanding way of isolating the partonic 
part accessible to pQCD from the non-perturbative parts. The basic ingredients are: (a) 
the hadron DA 0, which encodes the non-perturbative information regarding the momen- 
tum distribution of the constituent "near" on-shell valence partons collinear to the hadron 
and also features of the QCD vacuum structure as expressed through the quark condensates 
0,0,0]' an d (b) a scattering kernel Th, describing the hard scattering of "far" off-shell 
valence partons. The overall amplitude of the exclusive process is then given by the convo- 
lution, (p® Th ®4>. However, the application of pQCD to exclusive processes at intermediate 



momentum transfers or phenomenologically accessible energies (e.g ., at CEBAF) 



78 



79|,|80, 



81 



ras been 



82 



83 



84 



85|. It is 



the subject of severe controversies and criticisms [69l . 
widely anticipated that non-perturbative effects arising from soft gluon exchanges or from 
endpoint contributions to phenomenologically acceptable wave functions (or DAs) dominate 
and may severely preclude the predictability of pQCD. Hence, in this paper we use a modi- 
fied "resummed" pQCD formalism (as proposed in Refs. 22], |23J]) which is believed to largely 
enhance the predictability of pQCD in a self-consistent way at intermediate energies. The 
central issue here is the inclusion of transverse momentum kx dependence that necessitates 
the inclusion of a Sudakov suppression factor. This is to organize the large double loga- 
rithms of the type a s In 2 kr, arising at all orders due to the overlap of soft and collinear 
contributions of radiative gluon loop corrections. Such a resummation effectively suppresses 
the non-perturbative contributions at large energies. However, this may not still be effective 
enough when talking of Q 2 down to a few GeV 2 . These facts are in agreement with some 



of the recent findings, reported in Refs. 62 
context of B systems [86 ] . 



63 



64j for the pion form factor and also in the 



In this paper, we emphasize the importance of including two distinct contributions to ex- 
clusive quantities for obtaining good agreement with experimental data at low and moderate 
energies: firstly, the non-factorizable soft contributions which are not calculable within the 
perturbative framework and secondly, the power suppressed corrections from non-leading 
twist structures (twist-3) determining the preasymptotic behavior. In other words, the elec- 
tromagnetic form factor Fm{Q 2 ) for a charged meson M should be written as 



Fm(Q ) = FT(Q 



-ihard 

M 



(Q 2 ) 



(i) 



where F^f rA {Q 2 ) is the factorizable part computable in pQCD and F|° ft (Q 2 ) is the non- 



factorizable soft part. The soft contributions to the form factors can be calculated using 
phenomenological quark models either incorporating transverse structure (momentum) de- 
pendence of the hadron wave functions (see e .g. , Refs. 



rules via Local Duality (see e.g., Refs. 



48 



49 



50 



51 



61 



6l78 



80 



81 



87] ) or from QCD sum 



88j.) In the present paper, we follow 



the latter approach. We also focus on the presence of non-perturbative enhancements arising 
from kinematic endpoint regions of the scattering kernel which tend to invalidate collinear 
factorization. For our calculations, we use model twist- 2 and twist-3 light-cone wave func- 
tions incorporating transverse degrees of freedom, where the collinear DAs are derived from 



QCD sum rules [34], |35|, |36j, |37j, 138]. Naively, the twist-3 contributions to the form factor 
are expected to be small compared to leading (twist-2) contributions as they have a relative 
1/Q 2 suppression. On the contrary, the existing literature, either using model or asymptotic 



DAs 



25 



27 



28 



57 



641 ]. shows large twist-3 corrections to the pion form factor which 



even overshoot the twist-2 contributions in a wide range of low and intermediate energies. 
This is also confirmed in our analysis and is in fact more enhanced for the kaon form factor. 
To this end, our analysis shows good agreement with the existing pion data and in addition 
we prove the consistency of our results by adopting a scheme of analytization of the running 



strong coupling 



61 



J, |90J that removes the explicit Landau singularity at Q 2 = Aq CD 
by a minimum power correction in the UV regime. 

The paper is organized as follows: In Section [Til we briefly discuss the idea of factorization 
and review the basic definitions of the twist-2 and twist-3 pseudo-scalar meson DAs and 
their renormalization evolutions. Section II I II deals with the theoretical framework involved 
in calculating the space-like electromagnetic form factor. Here, we recall the predictions of 
classic asymptotic QCD for large Q 2 — >oc and how one needs to modify pQCD with collinear 
as well as "fcy" factorization schemes including Sudakov effects at intermediate energies. In 
Section \IV\ we provide the details of our numerical results for the pion form factor and 
compare it with the available experimental data. We also give a preliminary prediction for 
the kaon form factor, despite the lack of available experimental data for comparison in the 
desired phenomenological regime. Finally, Section |V] contains our summary and conclusions. 
The appendices contain a compendium of relevant formulae used in our analysis. 
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II. FACTORIZATION AND DISTRIBUTION AMPLITUDES 



The parton model of describing exclusive processes in QCD inherently rest on the so- 



called frozen approximation 



m 



71 



72j]). At high energies, exclusive scattering amplitudes 
are dominated by hadronic Fock states with essentially valence quark configurations (qq in 



mesons). While the relative velocities of the participating hadrons are located closely to the 
null-plane, the internal hadron "quantum-bindings" processes are highly time-dilated with 
respect to the exclusive reaction time scales in the rest frames of the remaining hadrons. 
This effectively freezes the hadronic internal degrees of freedom as seen by the other hadrons. 
This incoherence between the long-distance intra-hadronic binding processes and the short- 
distance inter-hadronic scattering reaction is the very motivation for the idea of factorization. 
Thus, the hadrons may be considered to be consisting of definite valence quark states denoted 
by a DA of leading twist 0. The collinear factorization formula is then used to express 
exclusive quantities like the form factors as a convolution using the DAs: 



F M {Q 2 ) = / dxdy(p in (x,iil)T H (x,y,Q 2 ,fil,^) 4> ont (y,fil) + ■ 
Jo 



(2) 



where Q 2 = — 2P in -P out . Here, P m and P out are, respectively, the ingoing and outgoing hadron 
momenta, x and y are the longitudinal momentum fractions of the nearly on-shell valence 
quarks, /^r is the renormalization scale and fiF is the factorization scale which is defined 
as the scale below which the QCD dynamics are non-perturbative and remain implicitly 
encoded within the DAs, while the dynamics above are perturbative and must be retained 
in the hard kernel Tjj. The ellipses in the above equation represent contributions from 
higher order Fock states and sub-leading twists which are all suppressed by inverse powers 
of Q 2 . In addition, they also include the non-factorizable soft contributions. Formally, the 



definition of the leading twist-2 DA for pseudo-sca 



34 



,|35|, 



70 



ar mesons (e.g., 7r ) can be given in a 



71 



72j in terms of matrix elements of 



process- and frame-independent manner 
non-local light-ray operator along a certain light-like direction z^ (z 2 = 0): 

(0\u{z)[z,-z]w B d(-z)\ir-{P)) = iP li [ dxe^ sp ^(x, t il)^ = 2x-l, (3) 

Jo 

with the path-ordering (V) Wilson line in terms of the gluon "connection" along the straight 
line joining z and —z along the null-plane given by 



z, -z 



V 



ig s / dy^A^y) 



(4) 



6 



where P 2 = m 2 and p^ is a light-like vector, 

The local limit z — > gives the normalization condition at an arbitrary scale fi, 

(p2-Tr(x, U 2 ) dx = 7^ (6) 

with the pion decay constant, f n ~ 131 MeV defined by 

(0|u(0) 7 ^(0)|7r-(P)> = i/ w P M . (7) 

The leading twist-2 DA 4>2;^{x, /x 2 ) can be expressed as a conformal series expansion over 
Gegenbauer polynomials C 2 3 ' 2 



"2k 



02 ;7r (x,/i 2 ) = -^x(l-x)(l + ^aL(// 2 )^ 2 (e)) , (8) 



n=l 



where 



r(as)/„\ j. /_ ..2 __\ ^/t 



S (z)=^2 ;7 r(z,M ^oo) = -^x(l-x) (9) 

is generally referred to as the asymptotic DA. The Gegenbauer moments a\ n represent the 
non-perturbative inputs encoding the long-distance dynamics and may be obtained e.g., via 
lattice QCD calculations or QCD sum rules. The renormalization group (RG) equation for 



< / ) 2;7r(a ; , H 2 ) is known as the Efremov-Radyushkin-Brodsky-Lepage (ER-BL) equation [70, 



71 



d r 1 



V -J^<P2-A X 1 V ) = J d V V ( X : V> )) <l>2;n(y, fl ) (10) 

with the integral kernel V(x, y\ a s ) to leading order in a s given by 



' " ' f l + -^—] 6{x -y) + -(l + -J— ) 6(y - x) 



V (x,y;a s ) = C F ^- .,, 

27r [ 1 — y \ x — y J y \ y — x 

where the "+" distribution is defined as 

[V(x,y;a s )) + = V(x,y;a s )-5(x-y) dtV(t,y;a s ) . (12) 

Jo 

Solving the above set of equations yields the multiplicative renormalization formula for 
moments a7 n to leading-logarithmic accuracy, 

a n (n 2 ) = L^l^aM) , (13) 



where L = a s (/i 2 )/a s (/io) and (3q = (11N C — 2Nf)/12, while the lowest order anomalous 
dimensions are given by 

^ =C F (i,{n + 2)+i,{l) 3 ' 



4 2(n + l)(n + 2)J ^ 

with the logarithmic derivative of the Gamma function ip(z) = T'(z)/T(z). Note that for the 
pion all odd moments a^ =135 ... vanish due to isospin symmetry. In contrast, the kaon DA 
have non-zero values for the odd moments signifying flavor-SU(3) violation effects. Hence, 
the twist-2 DA for the kaon (e.g., K~) is given by the expansion 

3/* 



V2Nc 



n=l 



f. 



K 



(j)2;K(x,fi 2 )dx-- 
l\j lPi c 



where the kaon decay constant fx ~ ^-22f n [24j is defined by 



(O |^(0) 7m 753(0)| K-(P)) = if K P, 



(15) 
(16) 

(17) 



The Gegenbauer moments being multiplicatively renormalizable with growing anomalous 
dimensions, for sufficiently large renormalization scale a finite number of moments are rele- 
vant, albeit the fact that the higher order moments have large uncertainties in their present 
determination. Hence, in all practical calculations, the series expansion of the DAs are trun- 
cated only to the first few moments. In this paper, we have adopted a model for the twist-2 
DAs in truncating up to the second moment, as was done in Refs. 35|, l37 |. 

For the charged pseudo-scalar mesons at the twist-3 level, there are two 2-particle DAs 
and one 3-particle DA. Here, we only give the formal definitions of the 2-particle DAs that 



we need in our analysis. For the charged pion (e.g., 7r ), they are defined by [34( : 



(0\u(z)i l5 d(-z)\n-(P))=fi n / dxe^^x,^) 



(0\u(z)a a ^f 5 d(-z)\ir-(P)) = -- f i 7r (P a z l3 -P l3 z a ) / dxe^</>^(x,fi 2 

Jo 



;is) 



with fi n = rri^/{m u + m d ) and similarly for the charged kaon (e.g., K ) [37J : 



Q\u(z)i7ss{-z)\K-{P)) = n K / dxe^(j>l K {x^ 2 ) 



(0\u(z)a a p l5 s(-z)\K-(P)) = -- f i K (P a z p -P p z a ) / dx e^<f>% K {x, fi 2 ) 

<J Jo 



(19) 
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with \lk = m 2 K /(m u + m s ). Note that the gauge-link factors (Wilson line (j3J)) in the matrix 
elements are to be implicitly understood. The twist-3 DAs have the following asymptotic 
forms: 



%S\x) = ^=x(l-x) ; M = «±,K± (20) 



with the normalization condition, 



o 



1 = ^ 



For our analysis, we use the 2-particle twist-3 DAs from Refs. [35|, |37|| defined at the scale 
fi — 1 GeV. As a matter of book-keeping, we explicitly provide the relevant formulae for the 
charged pion and kaon DAs in Appendix [A] 

III. SPACE-LIKE ELECTROMAGNETIC FORM FACTOR 

The electromagnetic form-factor is considered as an important observable for studying 
the onset of the perturbative regime in exclusive processes. For large Q 2 , the asymptotic 
scaling behavior Fm(Q 2 ) ~ ^/Q 2 follows from the well-known dimensional "quark counting" 



while for smal 
model 



42 



43 



Q 2 , the behavior is well described by the Vector Meson Dominance (VMD) 



44j and given by 

Fm{Q 2 ) « 1 + g2 1 //x 2 ^ Q 2 « ^vdm , (22) 

where //vdm ~ 750 MeV is a reasonable cut-off mass scale, showing no obvious trace of pQCD 
scaling behavior where there exists no high energy cut-off. Hence, a thorough understanding 
of this transition behavior (from nonperturbative to perturbative) is of crucial importance 
in QCD for understanding the very nature of strong interaction and in providing a vivid 
picture of the underlying quark-gluon substructure of the mesons. 

For a charged meson M (e.g., 7r ± , K ), the form factor is specified by the following matrix 
element: 

(P' + P),F M (Q 2 ) = (M(P') 1^(0)1 M{P)) ; J, = J2 e WWf , (23) 

/ 

where J M is the electromagnetic current with quark qf of flavor / and charge e/. In this 
paper, we shall only consider space-like momentum transfers i.e., q 2 = (P' — P) 2 = —Q 2 . 



Neglecting the meson masses, we consider the "brick wall" frame where the incoming particle 
with 4-momentum P in the z direction recoils with 4-momentum P' in the — z direction after 
interacting with the hard photon "wall". In the light-cone formalism, P = (Q/y/2,0,0 T ) 
and P' = (0,Q/y/2, Or). 



A. Hard contribution in pQCD 



The hard contributions to the form factor are calculated using the collinear factorization 
formula Eq. (j2J), where the hard scattering kernel T# at the scale // = /i p = //r is given to 
the leading order in a s by 

~2 1 1 1 



T H (x,y,Q ,//') = l<6-KC F a s (jJi 



(24) 



_3xyQ 2 3 (1 - x)(l - y)Q 2 _ 
where in QCD the value of the Casimir operator in the fundamental representation of 

SU(3) is Cp = (N 2 — 1)/2N C = 4/3. The factorization formula then yields the classic pQCD 

expression for the meson form factor at /i 2 = Q 2 : 

16nC F a s (Q 2 ) f 1 J <f)2-M(x,Q 2 ) ° 



T7'hard/' / o2\ 



Q 2 



dx- 



x 



(25) 



Note that using the asymptotic twist- 2 DA (p^f^x), one obtains the familiar 1/Q 2 scaling 



behavior for Q 2 — ► oo 



T7'hard/ / o2 N \ 

r a 



M 



(Q 



l7ia s (Q 2 )f 2 M 

Q 2 



(26) 



The principal motivation of the modified "resummed" pQCD is the elimination of large 
logarithms in the hard kernel that arise from radiative gluon loop corrections. One way of 
doing this is by the introduction of intrinsic transverse momenta dependence of the con- 
stituent partons, giving rise to a Sudakov suppression due to certain partial resummation 
of transverse terms, as mentioned earlier in the introduction. Including the transverse mo- 
menta of the two valence quarks within the meson, the tree-level hard kernel Tjj in the 

momentum-space is written as 

2 2 lQnC F a s (ii 2 )xQ 2 



T H (x,y,Q ,k 1T ,k 2T ,/i 



(xQ 2 + k 2 1T )(xyQ 2 + (k 1T -k 2T ) 



(27) 



where the transverse momentum dependence now sets the factorization scale. Then the 
modified factorization formula in the transverse impact parameter representation is given by 

d 2 h d 2 b 2 



7 hard 
M 



(Q 



dxdy 



J (27r) 2 (2n) 2 
xV 2; m{x, h, P, y) T H {x, y, Q, b u b 2 , n) V 2;M {y, h, P' , fj) 



(28) 



10 



where the modified DA V 2; M(xi, h, P%, A*) absorbs the large infrared logarithms into the 



Sudakov exponent S{ 



221 ] (including also the evolution of the DA from the factorization 



scale 1/bi to the scale p): 



V 2;M {xi, h, Pi ^ Q, /■«) = exp [Si(XQ, b h /i)] V 2 - >M {xi, h, 



s(AQ, 1/6, 



Si{XQ, b u y) = s( Xi Q , 1/6;) + s((l - a*) Q, l/h) + 2 / ^ lq (a s (jl 2 )) 

Jl/bi H 



(29) 



h/bi 

where l/b\, l/b 2 set the factorization scales in the transverse impact configuration. In the 
above equations, the quark anomalous dimension is given by 7 g (a s ) = —a s /ir and the "cusp" 
anomalous dimensions A and B, to one-loop accuracy are given by 



+ 



7T 



67 ^ 



10 



2a s (/i 5 



3 7T 



In 



JV C Nf + -3 ln[ 

27 9 J 27 i 3 M0 V 2 



,7b 



a s (/i 



2 



7T 



(30) 



where the MS running coupling to two-loop accuracy in standard perturbation theory is 
given by 

« s (/i 2 )_ 1 Aln(ln(/. 2 /A 2 QCD )) (0>i) 



Aln(ln(/. 2 /A 2 QCD )) 



7r /3 ln(/x 2 /AQ CD ) Mo m VF /iv QCDy 
with #j = (llA / c -2AT / )/12 =9/4 and (51iV c - 19JV»/24 =4 for N c =N f = 3. Note that 
the above modified factorization calls for introducing a scale hierarchy XQ > 1/bi > Aqcd 
(where X = Xi, (1 — xi), X\ = x and x 2 = y) to separate the distinct contributions from the 
perturbative and non-perturbative kinematic regions without the possibility of a "double 
counting". Note that there exist other schemes of defining the running coupling involving 
power corrections, restoring the explicit Landau singularity and the analyticity at Q 2 = 
(see e.g., Refs. 



89 



90] and also Section [TV] for details.) 



At low momentum transfers, the modified infrared free DAs are often approximated with 
constituent quark masses which are different from the actual masses of the current quarks 
and usually chosen close to the intrinsic transverse scale Aqcd of the hadron structure, i.e., 
between 200-500 MeV. These quark masses which effectively parametrize the QCD vacuum 
effects are also used to suppress possible endpoint effects. Hence, we have 



V 2 -M{xi, h, l/h) ^ V 2; M{xi, h, 1/bi, 



(32) 
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which could be expressed in terms of the full momentum-space light-cone wave function 
^2 ; m (which also includes the transverse momentum distribution of the constituent bound 
state partons): 



V 2;M {xi,K l/bi,M q ) 



d 2 k lT 

- • , 1 //, ,i 167T 3 



^ 2 . M (x i ,k iT ,l/b i ,M q ) . 



(33) 



/k2 T <(l/60 2 

To model the intrinsic transverse momentum dependence of the meson wave functions, we 
use the Brodsky-Huang-Lepage (BHL) gaussian prescription 



81 



^ 2 ;M(xi, ki T , M q ) = $ 2 ;M(xi, l/6») Ufa, k iT , M q ) 



(34) 



with 



$2;M(^, 1/bi) = A 2 -M 02;M(^i, V&O , 



16tt 2 /£ m r A 2 



( 1 



cxp 



2:M 



(14 + .M 2 .) 



(35) 
(36) 



assuming equal masses of the two constituent quarks within the meson. The parameters 
A 2 -,m, P2;M and 7W g are fixed using phenomeno logical constraints. The above integration 
then yields the full modified wave function in the impact representation: 



V 2 -M(xi, k, M q ) = A 2 - M 02;m(^, 1/bi) exp 

bfxiil - Xi) 



Xi (1 Xi ) 



xexp 



(37) 



Including the RG evolution equation for the hard kernel: 



f H (x,y,Q,bi,b 2 ,ii) = exp 



-4 f^lMsif)) 



f H (x,y,Q,bx,b 2 ,t) 



where 



t = max(y/xy Q, l/b\, l/b 2 ) 



(38) 



(39) 



one arrives at the "double-b" factorization formula for the meson form factor at the twist- 2 
level (5^] : 

p 1 ^00 

Fm =2) (Q 2 ) = 16vrQ 2 C F / xdxdy / Pa^fo &i, V&i, 



xV 2 - M (y, l/b 2 , M q ) H(x, y, Q, 61, 6 2 ) exp [S(x, y, b h b 2 , Q)} 



(40) 
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with 



H(x,y,Q,b 1 ,b 2 ) =K (^Qb 2 ) [9{b x - b 2 )K (^Qh)I (V^Qb 2 
+ 6(b 2 - b 1 )K ( v ^Qb 2 )I ( v ^Qb 1 ) 



(41) 



Kq and Jo are modified Bessel functions and the full Sudakov exponent is given by 

2 



S(x,y,bx,b 2 ,Q) = 



i=l 



s( Xl Q, l/bi) +s((l- Xi ) Q, 1/h) + 2 / ^t lq (a s (ji 2 )) 

1/6, P 



(42) 



For completeness, the expression for the full Sudakov factor s(XQ,l/bj), up to next-to- 
leading logarithm accuracy is given in the Appendix [Bl The expression slightly differs from 
the result given in Ref. 9lj], but numerically this difference is insignificant at our working 
accuracy. Note that this difference was first observed in Ref. 6o| . 

To include the sub-leading twist-3 corrections to the form factor, the hard scattering 
kernel gets slightly modified as compared to the twist-2 case which turns out to be [51 

647tCf a s (/i 2 ) x 



Th 3) (£ 5 2/ ; <3 2 ; ki T ,k 2 T,/i 2 ) — -, - , , > , , , 

(xQ 2 + kf T )(xy Q 2 + (kiT - k 2T ) 2 J 

Applying the momentum projection operator 



(43) 



92 



24 



V ■ u — v 3;M . ■ D u^3;Af O 



(44) 



a/3 



on the bilocal matrix element with quark flavors fx and f 2 (f\ 2 = u, d, s) 



(0\q f M<lf 2 (-z)\M(P)} a p = i 



dx / d 2 k T e^ zp) 



(45) 



ai3 



where ^%'. M {x, k T , 1/6, M q ) = d^. M (x,k. T ,l/b,M q )/dx, n = (1,0, T ) is the unit vector 
in the "+" direction, n = (0, 1,0t) is the unit vector in the "-" direction, £ = 2x — 1 and 
/j,m — m2 M/{ m q f + m q f ); one obtains the final formula for hard meson form factor up to 
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twist-3 corrections given by |62l. \Q4 \ 

Fir d {Q 2 ) = Ftr 2 \Q 2 ) + Ftr z \Q 2 ) 



p 1 /*oo 

32ttQ 2 Cf I dxdy / bxdbib 2 db 2 a s (t) 
Jo Jo 



x 



x ~ 



V 2 - M (x, b u l/b u M q ) V 2 , M {y, 6 2 , l/6 2 , M q ) 



+ ^7 {xVi; M (xM,Vb l ,M q )Vg. M (yMA/b2,M q ) 
Hm& h, l/h,M q ) P 3 ° M (y, b 2 , 1/62, A4 9 ) 



+ 



1 ~ 



+ 2 P 3-Af (x, 61, V&i, -H) ^ M (y, 63, l/6 2 , A4 9 



x 



H{x,y,Q,b 1 ,b 2 ) 

Y\_St{xi) S t (x i )exp[-S(x,y,b 1 ,b 2 ,Q)] ; x; = 1 - X; . (46) 



8=1 



Here, we have assumed a similar gaussian ansatz in the transverse momentum distribution 
of the modified twist-3 wave functions: 



Vl M {x u bi, 1/6;, M q ) = A P 3 . M (f>l jM (xi, exp 

b 2 Xj(l - x, 



X{ (1 X;) 



xexp 



V£ iM (xi, k, 1/6;, Mg) = Al M (j)% M {x h 1/6;) exp 

6, 2 X;(1 - X;) 



Xi (1 X 7 ; ) 



xexp 



(47) 



The hard kernel if and the Sudakov exponent S are given by Eq. (1411) and Eq. (1421) . re- 
spectively. The above formula is used to evaluate the pion and kaon hard form factors 
using the twist-2 DAs Eq. (JBJ) and Eq. ffT5]) . respectively, and twist-3 DAs provided in Ap- 
pendix [A] The S t {xj) are jet functions, defined as eikonalized matrix elements of quark 
fields attached by a Wilson line, arising from another kinematic resummation scheme called 
the threshold resummation, as introduced in Ref. 9J, [95|. The modified treatment of the 



collinear factorization prescription works reasonably well for the twist-2 case but for the 
twist-3 case, the Sudakov suppression factor may still not be effective enough in shielding 
the non-perturbative enhancements due to endpoint singularities. These are kinematic sin- 
gularities of the scattering amplitude when the longitudinal momentum fraction x of the 
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valence partons (quarks) go to 0, 1. Therefore, in addition there is a need to sum up the 
collinear double logarithms of the type a s In 2 x to all orders, which are then collected into 
these jet functions. The exact form of S t {xi) involves a one parameter integration, but for 



the sake of numerica 
proposed in Ref. [94 . 



calculations it is convenient to take the simple parameterization, as 
2 1+2 T(3/2 + c) 



95| : 



S t (xi 



\Xi ( 1 %>[ 



(48) 



VtFT(1 + c) 

where the parameter c ~ 0.3 for light pseudo-scalar mesons like the pion and kaon. The jet 
functions vanish at the endpoints and modify the endpoint behavior of the DAs, providing 
enough suppression to damp the artificial effect of endpoint singularities. 



B. Soft contributions via Local Duality 



The perturbative predictions for the pion form factor are known to be relatively sma 



48 



49, 



50 



69 



78 



83 



1 for 

as 



phenomenological low momentum transfers (Q 2 < 10 GeV 2 
also evident from our analysis in the next section. Clearly, there is the need for including non- 
factorizable soft contributions to explain the experimental data. The factorization ansatz 
Eq. (J2J) holds for large momentum transfers under the assumption that only the contributions 
from valence parton states dominate. This approximation no longer holds true at small 
momenta when contributions from higher Fock states with more than valence partons become 
significant. In addition, there could be non-perturbative enhancements from the so-called 
Feynman mechanism, which corresponds to selecting a hadronic configuration in which one 
of the valence parton carries almost the entire hadron momenta. Unfortunately, due to 
the complexity of soft QCD processes, there are no unambiguous ways to calculate these 
contributions analytically using the parton picture and Feynman diagrammatics, other than 
using theoretical models for the DAs. In this paper, we follow the Local Duality (LD) 
approach from QCD sum rules as in Ref. [61J where the same problem is addressed without 
a direct reference to DAs. In this section, we simply use the result for the soft form factor 
derived in the LD approach: 

l + 6s (QW 



rrrsoft / /^\2 \ 

r a. 



M 



Fl?{Q 2 ) 



(i + 4s (Q 2 )/Q 2 ) 3/2 ' 



(49) 
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The duality interval so, encodes the non-perturbative information about higher excited states 
and continuum contributions and is given by 

*o(Q 2 ) = 4vrW (1 + . (50) 

7T 

Expanding in inverse powers of Q gives F^f t {Q 2 ) ~ 1/Q 4 for large Q 2 and is thus expected 
to be sub-leading compared to the leading perturbative contribution from Eq. ( l26l) . Never- 
theless, at low and moderate momentum transfers the soft contributions turn out to be so 
significant in obtaining a good agreement with the experimental data. This fact is clearly 
revealed in our analysis in the next section. 

Next we add together the hard and the soft contribution to obtain the total contributions 
to the electromagnetic form factor Fm{Q 2 )- Here, it is necessary to ensure that the respective 
contributions lie within their domains of validity to minimize the possibility of a double 
counting. This technique, as introduced in Ref. 61[, employs gauge invariance that protects 
the value -Fm(O) = 1, through the vector Ward identity relating a 3-point Green function 
to a 2-point Green function at zero momentum transfer, i.e., F^f 1 (Q 2 = 0) = 1 . This implies 
that F\f vd (Q 2 = 0) = 0. A "smooth" transition from the hard to the soft behavior is then 
ensured by a matching ansatz from the large Q 2 behavior (arising from F^ (Q 2 )) to the 
low Q 2 behavior (arising from F|° ft (Q 2 )). This can be done by introducing a mass scale Mq 
which in the LD approach should be identified with the threshold Mq = 2sq. The twist-2 
part of the hard form factor F^j~ 2 \Q 2 ) is then modified following Ref. 61] 

2^Tv) F ^- < 51 > 

However, for the twist-3 case, the "matching function" &(z) = 1/(1 + z) 2 , with z = Q 2 /Mq 
is insufficient to ensure the Ward indentity at Q 2 = 0. To correct for the singular (~ 1/Q 4 ) 
behavior, we make a similar modification of the twist-3 part via the replacement 

Fr'w 2 ) = Ftr"\Q^ - p Or 3) (Q 2 )j0Qi (52) 

with the choice of the matching function &(z) = 1/(1 + z 2 ) 2 . This yields the Ward identity 



corrected twist-3 part: 

As 2 {Q 2 ) + Q A 



Ft 3) (Q 2 ) - LJT tz ) 2 Ft 3) (Q 2 ) ■ (53) 
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Finally, we arrive at our expression for the total electromagnetic form factor for a charged 
meson M (n ± , K^), valid for all values of Q 2 and is given by 

Fm{q ] - 1 - ^T^rniow (54) 



2s (Q 2 ) + Q 2 / W ) + UW) + Q 4 1 Fm (Q 



where f£" 2) (Q 2 ) and f£" 3) (Q 2 ) are given from Eq. ( HOI) and Eq. (14=61) . respectively. 



IV. NUMERICAL RESULTS 



At first, we need to determine the pion and kaon gaussian parameters A 2 -m-, Aj.^, -Af.j^ 
and /5 2; m, f° r the twist- 2 and twist-3 light-cone wave functions, respectively For 

the pion, they are obtained from the two constraints: firstly, by virtue of the leptonic decay 
7i — > yUz/ M , we have the condition 



dx 



d 2 k T 

167T 3 



leading to 



and secondly, from 7r c 



* w (a?,k T , M u> , 



01 Ml. 



A v / dx<p n {x)exp 

Jo L x{l-x) 

77, we have the condition 



2 V ^2N~ C 



2VQ 



1 fn AT 

dxV v (x,k T = 0,M u , d ) ■ V " 



which implies 



exp 



x(l — x) 



V6 



(55) 



(56) 



(57) 



(58) 



/„ x(l-x) 

where we use the constituent quark mass A4 u ,d — 0.33 GeV for both the u and d valence 
quarks in the pion. In the case of the kaon, firstly from the leptonic decay K — > we 
have the constraint 



Jo 



-P. 



K 



M 2 M 



u.d 



1 — X 



K 



2^6 ' 



(59) 



As for the second constraint, no straightforward condition like Eq. (1571) could be obtained 
for the kaon. On the other hand, by virtue of SU(3) isospin symmetry, it is reasonable to 
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make an assumption that for the kaon the average transverse momentum-squared of the 
valence partons defined by 

JdxJ d 2 k T \\4\ \y K (x,k T ,M uAs )\ 2 



<4> 



K 



(60) 



JdxJ d 2 k T \^ K (x,k T , M u ,d,s)\ 
has about the same value as in the case of the pion. We have checked that for both the 

1 /2 

twist-2 and twist-3 pion wave functions (k|) J « 0.35 GeV. This yields our second condition 
for determining the wave function parameters: 

1 Jodx(f) 2 K (x)exv 

(0.35) 2 



Ml , Mia 
x 1—x 



-1 >l 1 



+ 



M 2 . 

u,d 

l-x 



(61) 



where M. s = 0.45 GeV is used as the constituent s-quark mass and the full light-cone kaon 
wave function is given by 



ty K (x,k T ,M u ,, 



16ti 2 P 2 k A k 



x(l — X 

with x being the longitudinal momentum 
we use typical "double 



K [x)exp 



k 2 T + M 2 | k 2 T + Ml 



(62) 



of QCD sum rules 



34 



lumped" type [7 



35 



x 1 — x 

raction of the s quark. For our numerical analysis 



7l|, [72| DAs 4> nt K{x, /-f 2 ), derived in the framework 



371 ] . Note that we have considered N c = 3 in the expressions for 



the DAs. In Fig. [Hand Fig. El we display the twist-2 and twist-3 light-cone wave functions for 
the pion and kaon, respectively, along with their corresponding asymptotic wave functions. 
Note that the plots exclude the normalization factors of and for the individual DAs 
to facilitate comparison with one another. All the DAs are defined at the scale no = 1 GeV. 



The twist-2 and twist-3 DA input parameters are taken from Table 3 of Ref. [37J which we 
again provide in Table [I] along with the rest of the input parameters for the wave functions. 
Note that for the kaon we have shown both the type of wave functions, i.e., with and without 
including the (^parity-breaking terms. 



In Refs. 



34 



35 



371 ]. the DAs were assumed to obey equations of motion (EOM) of on-shell 



quarks for which fiM = rn 2 M /(m q + m q ^ s ) ~ 1.7 GeV was used. This is not strictly correct, since 
the quarks are not exactly on-shell but instead confined within the hadrons. We therefore 
prefer using a "chiral-enhancement" parameter x 3M (lGeV) ~ fi w ~ fix, instead of \iu in 
both the DAs and also the expression for the hard form factor Eq. (T46l) . Its numerical value 
is fixed by fitting the total form factor Eq. fl54l) to the available "world-data" for the pion 

Note that in this fitting procedure only the asymptotic forms of 
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FIG. 1: Twist-2 light-cone wave functions for (a) the pion Vi-^ and (b) the kaon Vi-k (solid lines), 
along with the wave functions corresponding to the respective asymptotic DAs (dashed lines). The 
DAs are defined at the scale /io = 1 GeV. 

the twist-2 and twist-3 DAs (Eq. ([9]) & Eq. ( l20i) ) were used in the pion wave functions. 
The running behavior X3m(a0 is then later introduced while calculating the form factors 
whose RG behavior is assumed to be the same as that of hm (see, Appendix [A]). In other 
words, this amounts to the replacement n 2 M — > X3m(1/^i) X3m(1/&2) m Eq. f|46|) . Also, as 
the bulk of the world pion data is concentrated in the very low energy region where the 
usual running coupling rapidly diverges, we in addition use an analytic prescription for the 
QCD running coupling to prove our results. The analytic scheme was suggested originally 



in Ref. 



60 



6l\ for 



90] for calculating the pion form factor and further developed in Refs. 
NLO calculations. Here, the central idea is the removal of the explicit Landau singularity 
present in perturbation theory rendering the coupling constant IR stable and reducing the IR 
sensitivity of perturbatively calculated hadronic observables. The scheme is also known to 
display higher loop stability. Now, the usual two-loop running coupling Eq. ( 13TT) in standard 
pQCD can be approximately expressed via the Lambert W_i function 



7T 



Pi 



1 + 



1% 



Pie 



A 2 

ii QCD 



-, -1 



(63) 



The extension of the above formula in analytic perturbation theory is too complicated to 
be evaluated exactly and instead there is an alternate approximate expression in the MS 
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FIG. 2: 2-particle twist-3 light-cone wave functions for the pion (a) T%.^ & (b) V%. m and for the 
kaon (c) T%, K & (d) f%. K with G-parity even terms (solid lines), along with the wave functions 
corresponding to the respective asymptotic DAs (long dashed lines). For the kaon, the twist-3 wave 
functions including G-parity odd terms are also shown (dotted lines). The DAs are defined at the 
scale fj,Q = 1 GeV. 



scheme, as suggested in Ref. [89| : 



a 



an^pprox/ . .2 



71 



1 

A) 



l l 

/ 1 — exp(Z) 



(64) 
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7r± 


At no = 1 GeV 




At hq = 1 GeV 


units 






fnu,d 


5.6 ± 1.6 [371 


MeV 






m s 


137 ± 27 [37] 


MeV 




0.33 


M u ,d 


0.33 


GeV 






M s 


0.45 


GeV 


m n 


139 


m K 


493 


MeV 


aj 





af 


0.06 ±0.03 [37] 






0.25 ±0.15 [37] 


af 


0.25 ±0.15 [37] 




u 


131 


Jk 


1.221 /J [24] 


MeV 


/3tt 


0.0045 ± 0.0015 [37] 


hx 


0.0045 ± 0.0015 [37] 


GeV 2 


^3tt 


-1.5 ±0.7 [37] 


^>3K 


-1.2 ±0.7 [37] 




^3tt 





^3K 


1.6 ±0.4 [37] 





TABLE I: Various input hadronic parameters for twist-2 and twist-3 light-cone wave functions at 
Ho = 1 GeV. 

where A an in the analytic scheme is the analog of Aqcd in the usual perturbation theory and is 
chosen to be around 0.4 GeV for Nf = 3. We use this formula for the analytic coupling in our 
calculations. The simple one-parameter fitting of the total form factor to the experimental 
data gives the best fit values of 1.2 GeV and 1.4 GeV for the usual and analytic QCD 
coupling schemes, respectively. Here, we choose the average value X3M = l-3 GeV for both 
the schemes, and generously consider the resulting difference from the phenomenological 
value of 1.7 GeV to contribute to the theoretical error, i.e., ±0.4 GeV. Note that the fitting 
takes into account the individual error-bars of the data points. 

Finally, following Ref. the Sudakov suppression factor exp(— s(XQ, 1/b)) is set to 
unity for small transverse separation "6" between the valence quarks, i.e., whenever b< 
a/2/ (XQ). Also, to avoid probing into certain kinematic regions where exp(— S) may become 
greater than unity causing an enhancement instead of a suppression, exp (-S) is set to 1 for 
S < 0. 
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7r± 


G (even) 


K± 


G (even) 


G (even + odd) 


units 


4. 
A %, 
(/V)' 

(/3f;J 2 


1.69(1.66) as 
3.76 (3.59) as 
3.37(3.33) as 
0.76 (0.87) as 
0.62 (0.74) as 
0.81 (0.87) as 


A 2 ;K 

A a 
3'K 

(PIk) 2 
(/?sV) 2 


2.06 (2.07) as 
4.40 (4.56) as 
4.16(4.14) as 
0.78 (0.89) as 
0.70(0.79) as 
0.88 (0.89)^ 


2.06 (2.07) as 
4.35 (4.56) as 
4.06 (4.14) as 
0.78 (0.89) as 
0.65 (0.79) as 
0.84 (0.89),, 


GeV" 2 
GeV" 2 
GeV- 2 


xt 


1.3 ±0.4 


xf K 


1.3 ±0.4 


1.3 ±0.4 


GeV 



TABLE II: Various determined hadronic parameters for twist-2 and twist-3 light-cone wave func- 
tions at fiQ = 1 GeV. The numbers in the parentheses correspond to values for the asymptotic wave 
functions. 

A. The pion form factor 



Using the DAs in Appendix A, we evaluated the total electromagnetic form factor for 
the pion Eq. ( 1541) using both the usual two-loop QCD running coupling Eq. ( 131~T) and the 
analytical prescription Eq. (1641) . Fig. [3] shows our results for the total form factor, along 
with the experimental "world-data" [l|, [2j, for the pion. The plots correspond 

to Aqcd = 0.2 GeV and A an = 0.4 GeV, respectively. It appears that the full twist- 
3 calculations improve the agreement with experimental data to a much better extent at 
intermediate energies down to around 1-2 GeV 2 than for the twist-2 case. Note that in the 
usual perturbative scheme, as Q 2 — > the total form factor becomes very unpredictable 
and starts oscillating between large values although F n (Q 2 = 0) = 1, satisfying the Ward 
identity. This clearly signals the breakdown of perturbation theory at such small momentum 
transfers. 

To study the contributions of endpoint effects and to distinguish individual soft and hard 
contributions, it is more useful to study the variation of the scaled pion form factor Q 2 F n 
with Q 2 . In Fig. HJ we show the individual contributions of the twist-2 and twist-3 power 
correction to the scaled hard pion form factor over a wide range of momentum transfers for 
the usual QCD coupling. Clearly, the twist-3 contributions are seen to be significantly larger 
than the leading twist-2 counterparts at low momentum transfers, supporting the claims 
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FIG. 3: The total electromagnetic pion form factor at the twist-2 level (Soft+Twist-2), denoted by 
the dashed lines and at the twist-3 level (Soft+Twist-2+Twist-3), denoted by the solid lines, with 
(a) the usual QCD running coupling and (b_^ the^analytical QCD running coupling. The world 
pion data are taken from Refs. 



upline and _[b) the ; 

Q,y,y,y,Hy,Q,3 



made in 



25 



27 



28 



57j, 



62 



64j. In fact, it is interesting to see the endpoint enhancement in 
the twist-3 amplitudes much more explicitly if one rather considered only the collinear DAs 
to calculate the form factors in the usual perturbation theory, without considering the full 
transverse momentum dependence (e.g., the BHL ansatz) in the meson wave functions, as 



95l |. In other words, one simply makes the replacement Vm{x) 



originally done in Ref. 

4>m{x) i n calculating the hard form factor. The inclusion of the transverse momenta and 
constituent quark masses in the wave function provides a natural cut-off for the soft and 
endpoint enhancements. Similar behavior can also be observed in the analytic case, although 
we have not displayed the enhancement which is less drastic. These facts suggest that the 
modified collinear factorization scheme, including explicit transverse degrees of freedom 
with Sudakov suppression, which works well for the twist-2 case is not very effective at the 
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FIG. 4: Twist-2 (long dashed lines) and Twist-3 (solid lines) corrections to the scaled hard pion 
form factor with the usual QCD running coupling. Plots (a) & (b) are obtained using BHL ansatz, 
while (c) & (d) are obtained with Vm{x) 4>m(x). Also, the plots (a) & (c) do not include 
threshold resummation in the hard contributions, which are included in (b) & (d). The soft 
corrections (short dashed lines) are also shown. 

twist-3 level in shielding such artificial non-perturbative enhancements at low momenta. To 
improve this situation, especially for the results obtained in usual perturbative scheme, we 
use threshold resummation which along with Sudakov suppression provide large damping of 
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the endpoint effects in the twist-3 amplitude. The twist-2 part on the other hand remains 
mostly unaltered, if not slightly enhanced due to the threshold resummation, especially in 
the low-energy region. Note that in this respect the use of threshold resummation in the 
analytic scheme is somewhat redundant and has little effect on both the twist corrections. 
Finally, as expected, one observes that the twist-3 corrections fall off rapidly with increasing 
Q 2 and beyond a certain point fall below the twist-2 corrections. At asymptotically large 
momentum transfers, only the twist-2 contributions are expected to dominate. 

Our final results for the scaled pion form factor are summarized in Fig. [5] and Fig. [61 
We use both the usual and the analytic QCD running couplings and compare our results 
with the available experimental pion data with increasing error bars towards intermediate 
energies. The individual soft and hard contributions along with the total contribution are 
shown. Clearly, contrary to the earlier claims made in Ref. (ssj, the twist-2 hard form 
factor is far too small in the phenomenologically accessible region to explain the data. One 
must therefore look for other possibilities like non-perturbative higher twist effects and soft 
contributions. Interestingly, it is seen that the soft dynamics largely dominate the low-energy 
region below 10-16 GeV 2 but rapidly fall off in the asymptotic region. The contributions 
from the twist-3 corrections turn out to be significantly large in the moderate range of 
energies below 100 GeV 2 , but eventually the hard twist-2 contributions solely determine the 
asymptotic trend beyond Q 2 « 100-150 GeV 2 . 

As evident from the figures, the total scaled pion form factor (solid lines) up to twist-3 
corrections displays an obvious improvement of the overall agreement with experimental data 
compared to the twist-2 scaled form factor (dashed lines). To some extent, it is somewhat 
surprising to see that in combination with the soft contributions, the modified resummed 
pQCD with the usual QCD coupling could work so well as low as Q 2 ~ 0.25 GeV 2 , far lower 
than previously envisaged. To this end, we display the results in the analytic scheme to 
confirm our results. The analytical prescription is known to reduce the scheme and renor- 
malization scale dependence, largely increasing the stability of solutions 6l|]. Accordingly, 
there is some confidence in our displayed results. The results obtained in both schemes 
not only show a striking similarity even at sufficiently low momentum transfers, but also 
show a good agreement with the available pion form factor data. However, whether or not 
such an agreement is merely accidental is matter of debate. There may still be substantial 
sub- leading contributions e.g., from a full NLO calculation in the strong coupling including 
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FIG. 5: Scaled total electromagnetic form factor for the pion with (a) the usual QCD running 
coupling and (b) the analytical running coupling in the low and intermediate energy regime. The 
solid line represents the full twist-3 result (Soft+Twist-2+Twist-3), the long dashed lines represent 
the twist-2 result (Soft+Twist-2) and the soft corrections are indicated by the short dashed lines. 
The world pion data are taken from Refs. 
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FIG. 6: Scaled total electromagnetic form factor for the pion with (a) the usual QCD running 
coupling and (b) the analytical running coupling over a wider range of the intermediate energies. 
The solid line represents the full twist-3 result (Soft+Twist-2 +Twist-3), the long dashed lines 
represent the full twist-2 result (Soft+Twist-2) and the soft corrections are indicated by the short 
dashed lines. The world pion data are taken from Refs. 

subleading twists and intrinsic transverse momenta for the hard scattering kernel and the 
DAs, or from higher order Fock states and helicity components in the light-cone DAs. Note 
that a NLO calculation in the strong QCD coupling constant was done in Ref. 60(] where the 



corrections to the twist-2 pion form factor was found to be quite large. A full NLO calcula- 
tion for subleading twists, including also the transverse momentum dependence, is however, 
still missing. In Ref. 63], contributions of higher helicity states were found to lower the 
total pion form factor significantly. Hence, without systematically taking all of these effects 
in account, which is beyond the scope of the paper, no definitive statement can be made as 
to how well our results agree with the data. Moreover, the available data itself has very low 
statistics at intermediate energy and is plagued with large uncertainties. It is therefore diffi- 
cult to give a proper theoretical error estimate of our results, when the asymptotic formalism 
is itself largely unreliable in the region of our interest. What we have done in this paper is 
a combination of model and pQCD calculation. It may thus be worth using the estimates 
for the ranges over which the input parameters of the DAs, namely fi w (x3n), /W, u 3n and 
a\ vary (given in the Tables HI fc HT|) . in computing our theoretical error. In addition, we 
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FIG. 7: Theoretical la-error for the scaled total pion form factor due to the variation of the input 
parameters for the DAs (x3ir, iW, ^3tt and a\ ) and AQCD.an- The lines correspond to the mean 
values of the parameters using the usual (dotted line) and the analytic (dashed line) QCD coupling 

schemes, respectively. The error bars for the experimental data points n y, y, y, y, y, y, 3 are 

also shown. 



allow a variation of ±0.05 GeV for both Aqcd and A an . We used a Monte Carlo technique 
to generate a gaussian "la" spread of the scaled form factor for various Q 2 values about a 
central mean. The error estimate displayed in Fig. [7] shows our maximum theoretical error 
to be about 10 percent for the usual QCD coupling and somewhat less for the analytic 
coupling. 
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FIG. 8: Scaled total electromagnetic form factor for the kaon with the usual QCD running coupling 
(upper plot) and the analytical running coupling (lower plot) for intermediate energies. The solid 
lines represents the full twist-3 result (Soft+Twist-2 +Twist-3), the long dashed lines represent 
the twist- 2 result (Soft+Twist-2) and the soft corrections are indicated by the short dashed lines. 
For the full twist-3 case, both the results i.e., with and without including the G-parity-breaking 
terms in the DAs are shown. 
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B. The kaon form factor 



We conclude the section on the numerical analysis by displaying our predictions for the 
kaon form factor, applying the same techniques as for the case of the pion form factor. 
The twist-2 and twist-3 light-cone wave functions for the kaon has the general form in the 
transverse 6-space given by 

b 2 x(l — x) 



Vk(x, b, 1/6, M U: d :S ) = A K <p K (x, 1/6) exp 

'M 2 



4# 



xexp 



-P. 



K 



+ 



X 



X 



(66) 



where we used the twist-3 chiral enhancement parameter Y3K = 1-3 GeV and the experimen- 
tal estimate for the kaon decay constant fx ~ 1.22f n 24j. Here, we also take Aqcd — 0.2 
GeV and A an = 0.4 GeV for the respective running couplings in the MS scheme. The results 
are summarized in Fig. [8] for intermediate energies. The solid line represents the total scaled 
form factor in each case. Here, our results for the kaon form factor must be considered 
preliminary. Due to the complete absence of experimental data at intermediate energies we 
are unable to make any meaningful phenomenological comparison. The presently available 
kaon data has very poor statistics and has hardly been measured above 0.2 GeV 2 . Hence, 
we do not show the experimental data points in the form factor plots. With the availability 
of better quality data in future there could be plenty of room for further improvements, 
for instance, extension of the above results to include a full NLO calculation for subleading 
twists or higher helicity and Fock state contributions. 



V. DISCUSSION AND CONCLUSION 



For the past two decades the electromagnetic meson form factors have been the subject 
of intensive theoretical and experimental scrutiny and yet there is still to be an universally 
accepted framework for their description. Presently, reliable experimental data are available 
only for the pions which are entirely concentrated at very low energies with very poor 
statistics at intermediate energies. The low energy part of the data is best explained by 
the standard VMD model, showing no apparent trace of pQCD behavior, which is expected 
only at very high energies. Very many attempts have been made to predict the onset of the 
perturbative behavior for the pion form factor. The modified or resummed valence pQCD 
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with factorization appears to show some attractive features to enable pQCD calculations 
to be valid in a self-consistent way even at very moderate energies. Whether this is true, 
will only be confirmed when data with better statistics at higher momentum transfers will 
become available in the future. At the same time, the onset of the perturbative behavior 
being very slow, it is still unclear whether the leading order perturbative calculations could 
be expected to be precise even at the highest accessible energies. In this paper, with the 
help of (a) the double-humped type DAs and (b) the modified transverse (kr) factorization 
scheme, incorporating both Sudakov suppression and threshold resummation, we got rid 
of non-perturbative end-point enhancements. This enlarges the scope of applicability of 
resummed pQCD independent of the coupling scheme to a much wider range of intermediate 
energies, if not down to a few GeVs, as demonstrated in this paper. By a simple adjustment 
of only the chiral enhancement parameter x 37r , a good agreement with the experimental data 
was obtained. 

As for the scaled pion form factor, we found that the leading order pQCD contributions 
potentially undermine the agreement with the available low energy data and are only trust- 
worthy in the hard energy regime: even from a very conservative point of view, Q 2 should 
be bigger than 4 GeV 2 . At low momentum transfers, the non-perturbative contributions 
dominate, being larger than the hard (twist-2) contributions at least by a factor of two. 
In fact below 4 GeV 2 , 60-70% of the available data is already accounted for by the soft 
contributions. In addition, we also needed the twist-3 power corrections to explain the re- 
maining discrepancy. However, at larger energies (say, Q 2 > 50-100 GeV 2 ), both the soft and 
the twist-3 contributions rapidly fall off and eventually the twist-2 form factor dominates 
asymptotically. Similar conclusions, albeit being preliminary, are drawn for the kaon form 
factor although it seems that the onset of the perturbative behavior occurs at slightly larger 
momentum transfers than for the case of the pion. Of course, as we mentioned earlier, it 
still remains to be investigated about the nature of the contributions that may arise from a 
full systematic NLO calculation with subleading twists and intrinsic transverse momenta, or 
from the inclusion of higher helicity and Fock states. With the availability of better quality 
of data in future such analyses may be necessary to make definite conclusions. 

In summary, although the quality of present experimental form factor data does not 
allow a definitive conclusion, one can expect that the non-perturbative soft contributions 
and higher twist power corrections to the form factors play an important role at phe- 
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nomenologically accessible momentum transfers. Thus, more work is needed to be done on 
both the theoretical and experimental sides to obtain more conclusive results and push the 
frontiers of our knowledge on confinement dynamics through the study of higher order and 
non-perturbative contribution to exclusive processes. 
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APPENDIX A: 2-PARTICLE TWIST-3 DISTRIBUTION AMPLITUDES 



The twist-3 DAs are obtained by an expansion over conformal spins. At next-to-leading 
order, the 2-particle DAs </>| M and M (including meson- mass corrections that break chiral 
symmetry at 0(m s + m q ) in the SU(3) case while preservin g G -parity) are given in terms 
the Gegenbauer polynomials ClJ 2 and Cn 2 , respectively 34, 35] as 



{l + (stW^) - \ pW)) cl /2 (0 



"2\ ^/m „\ 1 1 , { k„ f,.2\ 1„ ,..2, / 2, 



■\FW X<yl ~ X )y + [5V3m(V ) - 2^3Af(/i )W 3M (P ) 



with 



fzM 1 



Pm 



m M 



M = 7T ± , K* 1 



(Al) 



(A2) 



where the non-perturbative parameter f 3 M and c^3M, respectively are defined by the following 
matrix elements of local twist-3 operators: 



(0 \q h a^ b g s G a pq h \M(P)) = if 3M {P a P^g u p - P a P„g^ - PpP^g ua + PpPuQa^) , 



(°k/i<VA75 9sG a \)qf 2 - {?>/7)id l3 q h a flX ^g s G aX q h \M(P)) 
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= 14 ^ 3M ^' a ^^ uj3M ' (A-'V 

where G a p is the gluon field tensor. The LO scale dependence of various twist-3 parameters 
are given by 
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where L = a s (p 2 )/a s (p 2 ), C F = (N 2 - 1)/2JV C , C A = N c and p « 1 GeV. 

We have also considered the 2-particle twist-3 kaon DAs K and <p 3K , as given in ^t\ , 
which not only include a complete set of meson-mass corrections but also G-parity-breaking 
terms of 0(m s — m q ): 



with 



K 



1 + 3p^(l + Qaf ) -9p*af + 
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Cl /2 (0 + (sO V3K + lhpl4 - 3p A af )c 2 1/2 (0 
KWa* - \p-a(\ Cl /2 {i) - 3r) 3K u 3K Cl /2 (0 + | (p* + p5) 
(l - 3af + 6af ) Ins + ^ (p+ - p5) (l + 3af + 6af ) ln(l - x) J 



(A5) 



+ 



3/, 



2v/2iV: 



3 
2 



3p?af - H p ^af ) Cl'\) + ^ 3K - \ rj 3K co 3K + \ p K + a^ of (0 
+ V3k\ 3K C 3 3 /2 (0 + ~ (p+ + P K ) (1 - 3af + 6af ) lnx 
+ 1 (Pi - P K ) (1 + 3af + 6a*) ln(l - *) } (— L 



/W 1 K (m s + m g ) 2 K m 2 -m 2 

?73JC = , P+ = — and p_ ■ 



Ik Pk 



m 



K 



m 



(A6) 



(A7) 



K 



Note that the expression for (p^. K (x, p 2 ) is normalized to unity with an extra factor of 
1/(1 — p + ), compared to that given in [371 ] . The non-perturbative parameters f 3 x, ^3K and 
\ 3 k are defined (e.g., K~) by 

\ua^ 5 g s G a p s\K~(P)) = if 3K {P a P^g v p - P a Pv9^ - PpP^ua + PpP„g afl ,) , 



(0 | it cr^^iDp, g s G aX ]s - (3/7) id pilaris 9 s G a \s\ K (P)) = —if 3K P a PpP»U3K 



uiDpa^xIsgsGax s - ua^ 5 g s G a \iD p s 



K-(P)) = -tfsKPaPpP^SK- 
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where in the chiral limit the renormalization group equations at LO give, 



But the strange quark being massive, there is operator mixing of the ones in Eq. (1A8[) with 
that of twist-2 operators and the resulting LO renormalization group equations give the 
following scale dependence of the various twist-3 parameters: 
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Finally, we present the various Gegenbauer polynomials used in the above formulae: 



C 3 1/2 (0 = ^(5£ 2 -3), 
C 4 1/2 (O = ^(35£ 4 -30£ 2 + 3), 

c 3/2 (0 = i, 

C? /2 (0 = 3£, 

^3 /2 (0 = ^(7e 2 -3). (All) 



APPENDIX B: CALCULATION OF THE SUDAKOV EXPONENT 

The full expression of the Sudakov suppression factor S(x, y, bi, b 2 , Q) is given by, 



S(x, y, 61, b 2 , Q) = s(xQ, 61) + s(yQ, b 2 ) + s((l - ar)Q, 61) + s((l - y)Q, 6 2 ) 
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In the above formulae, 
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